NASA Contractor Report 4113 
DOT /FAA /CT-87 1 30 


Three-Dimensional Trajectory 
Analyses of Two Drop Sizing 
Instruments: PMS OAP 
and PMS FSSP 


Hillyer G. Norment 


CONTRACT NAS3-24630 
FEBRUARY 1988 


(NASA -C 5-4! 13) 1EBEE-C2BENS1C N A I IfiAJiSC'ICRi Nbfc-1657 

ANALYSES CE 'I WO DEOf* SI2ING J N Slfc U S : IMS 

CA i AND EMS FSSE final Fepcct (A tmcspnet ic 
Science Associates) 65 p csCi QIC unclds 

ill/OJ 0 1 2 b7 S 7 



NASA Contractor Report 4113 
DOT /FAA/CT-87 1 30 


Three-Dimensional Trajectory 
Analyses of Two Drop Sizing 
Instruments: PMS OAP 
and PMS FSSP 


Hillyer G. Norment 
Atmospheric Science Associates 
Concord, Massachusetts 


Prepared for 

Lewis Research Center 

under Contract NAS3-24630 


IW1SA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
Information Division 


1988 



FORWARD 


The able assistance of Angela Quealy, Sverdrup Technology, Inc., 
Middleburg Heights, Ohio, with part of the computing and plotting is 
acknowledged. 


PRECEDING page bunk not filmed 



TABLE OF CONTENTS 


Page 

SUMMARY 1 

INTRODUCTION 2 

THE INSTRUMENTS AND AIRPLANE 8 

THE INSTRUMENTS 8 

THE AIRPLANE 12 

MOUNTED INSTRUMENTS 12 

CALCULATION METHODS 17 

FLOW CALCULATION 17 

WATER DROP TRAJECTORY CALCULATION 18 

CONCENTRATION FACTOR 21 

RESULTS 24 

FLIGHT CONDITIONS 24 

FLOW RESULTS 24 

TRAJECTORY RESULTS 26 

OAP 35 

FSSP 42 

Preferred Orientation 50 

CONCLUSIONS 51 

ISOLATED INSTRUMENTS 51 

WING-MOUNTED INSTRUMENTS 51 

PREFERRED ORIENTATION 52 

APPENDIX A - CONCENTRATION FACTOR AND SPEED RATIO DATA 53 

APPENDIX B - DRAG VECTOR ANGLES 56 

REFERENCES 59 



SUMMARY 


Flow-induced distortions of water drop flux and speed for two hydro- 
meteor measuring instruments, Particle Measuring Systems OAP and FSSP, 
are predicted by three-dimensional calculations. The instruments are 
studied in isolation and as mounted under the wing of a DeHavilland Twin 
Otter airplane. Several free stream air speeds, and angles of attack of 0° 
and 4° are studied for drop diameters ranging from 2 to 1000 ym. 

For the OAP in isolation and under the Twin Otter wing with the 
airplane at 0° angle of attack, distortions of practical consequence are 
not found. At 4° airplane angle of attack, 17% undermeasurement of both 
flux and speed is predicted for cloud-size droplets. 

The FSSP presents greater flow obstruction than the OAP, and in addi- 
tion, air and drops must traverse the measurement tube. As expected, we 
predict larger flow-induced effects under all circumstances than for the 
OAP. For the FSSP in isolation and mounted on the Twin Otter at 0° 
angle of attack, both speed and flux are predicted to be undermeasured by 
about 10% for cloud-size droplets. At 4° airplane angle of attack, 24% 
undermeasurement of both flux and speed is predicted for cloud-size droplets. 

For both wing mounted instruments we find that a large portion of the 
flow induced effects, approximately half, is caused by the instruments them- 
selves. This shows that the common practice of neglecting instrument flow 
effects in comparison with aircraft flow effects can produce results that 
are seriously in error. 

Preferred orientation (canting) angles of distorted water drops as 
they pass through the instruments are predicted to vary with drop size, 
angle of attack and free stream air speed in apparently complicated ways. 

This needs further study. 
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INTRODUCTION j 

NASA Lewis Research Center has instrumented a DeHavilland DHC6 Twin 
Otter airplane (Fig. 1) for use in icing research. Instrumentation 
(Fig. 2) on the airplane is designed to collect data critical to aircraft 
icing during traverses through icing clouds. All of these instruments 
require the sensing or collection of water drops or other hydrometeors, 
and as is well known (refs. 1,2), flow perturbations caused by the pas- 
sage of the airplane through the cloud can act to seriously distort the 
measurements. Since experimental determination of these effects would 
be prohibitively expensive, three-dimensional flow and hydrometeor tra- 
jectory calculation methods (refs. 1, 2, 3) are used for their prediction. 

The instruments of concern in this study, the PMS* Optical Array Probe 
(OAP) and PMS* Forward Scattering Spectrometer Probe (FSSP) are numbered 
6 and 7 in Figure 2, and are mounted under the Twin Otter wing as shown 
in Figures 2, 3 and 4. Results of similar studies of other instruments 
shown in Figure 2 are presented elsewhere (ref. 4). 

There are two main objectives of this study: 

1. Predict distortions of water drop flux and speed caused by 
flow about the isolated instruments. 

2. Predict distortions of v/ater drop flux and speed caused by 

flow about complete assemblies of the instruments, their 
mounts and the Twin Otter airplane. 

Results are reported for water drops in the diameter range 2 to 1000 ym 
and for orientations of 0° and 4° angles of attack to free stream flows 
at several speeds. Also, results are compared with similar studies done 
by the Canadian National Aeronautical Establishment for their Twin Otter 
which is similarly instrumented (refs. 5, 6). 


* Particle Measuring Systems, Inc., Boulder, Colorado. 
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1. Leigh Ice Detector 

2. Pressure Ice Rate and Accretion Sensor (PIRAM) 

3. Rosemount Ice Detector 

4. J-W Liquid Water Content Sensor 

5. CSIRO-King Liquid Water Content Sensor 

6. OAP Laser Spectrometer 

7. FSSP Laser Spectrometer 

8. Experiment Carrier Port (Rotating Multi Cylinders) 

9. Soot Slide Droplet Sampler 

Figure 2. Mount sites of icing instruments on the NASA LeRC Twin Otter. 

Instruments considered here are numbered 6 and 7. Results 
of similar studies of the other instruments are described by 
Shaw, et al . (ref. 4) . 
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Only effects of flow perturbations on drop trajectories are considered. 
Measurement and sampling errors, and biases caused by optical problems or 
other workings of the instruments are not addressed. 
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THE INSTRUMENTS AND AIRPLANE 


THE INSTRUMENTS 

Figures 5 and 6 show computer plots of detailed digital descriptions 
of the instruments' surfaces suitable for use by the flow calculation 
codes. Dimensional data also are given in the figures. 

The optical array probe (OAP) (refs. 7, 8) consists of a cylindrical 
canister with hemispherical ends, one of which supports a pair of probe 
arms. A laser beam passes between the probe arms as shown in Figure 5, 
and hydrometeors that pass through the beam are sensed and sized. Two- 
dimensional images of the hydrometeors are recorded via use of a linear 
array of sensors, and to assure that undistorted images are obtained, 
transit speeds of the hydrometeors through the beam must be known. Maxi- 
mum diameter range capability of the instrument is approximately 75-2000 pm. 

The forward scattering spectrometer probe (FSSP) (ref. 8) uses a 
canister with hemispherical ends that is identical in size and shape with 
the one used by the OAP. Support arms, which are considerably larger in 
diameter than those of the OAP, extend from one end of the canister, and 
these support a measurement tube through which the particles must pass 
(Fig. 6). A laser beam passes diametrically across the inside of the 
measurement tube; particles are sensed as they intersect this beam and 
are sized by use of Mie scattering theory. Maximum diameter range capability 
is approximately 2-100 ym. 

Structural details of portions of the instruments that are located 
well aft of the sampling volumes contribute little of significance to flow 
at and forward of the sampling volumes, so that an abbreviated version 
of the canister was used for studies of the instruments in isolation. 

Figure 7 shows computer plots of the instruments with the abbreviated 
canister as used for the isolated instruments studies. 
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Figure 5. Detailed digital description of the PMS OAp suitable for use by the three- 

dimensional flow code. Dimensional data, labeling, and laser beam path lines 
are added to a computer plot. Probe arm diameter is 1", and laser beam 
diameter is taken to be approximately 0.7". 



MEASUREMENT TUBE 
SUPPORT ARM 



Figure 6. Detailed digital description of the PMS FSSP suitable 

for use by the three-dimensional flow code. Dimensional 
data and labeling are added to a computer plot. Inside 
diameter of the measurement tube is 1 V 16 ". Support 
arm diameter is 1 7 / 8 ". 
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b. PMS FSSP 


Figure 7. Digital descriptions of the PMS OAP and FSSP used for studies 
of the instruments in isolation. Aft ends of the cannisters 
are abbreviated; otherwise the paneling is the same as shown 
in Figs. 5 and 6„ 
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THE AIRPLANE 


A detailed digital description of the Twin Otter airplane has been 
prepared that includes all major components (Fig. 8). For this study this 
much detail is not needed, so that an abbreviated version of the paneling 
was used as shown in Figure 9. Since neglect of the engine nacelle might 
be questionable, concentration factor (defined below) calculations were 
done with the nacelle included for water drops to the OAP, which instru- 
ment is closest to the nacelle. Comparisons of these results with results 
obtained without the nacelle show unambiguously that the nacelle does not 
significantly influence the results. 

The wing airfoil section, NACA 0016, is spanwise uniform with a con- 
stant 78" chord. The wing has no geometric twist. To mount the wing on the 
fuselage the following operations are required: 

1. Rotation by 2.5° about a spanwise axis, leading edge upward. 

2. Dihedral rotation by 3° about a chordwise axis at a distance 
of 35.15" from the symmetry plane, tip upward relative to root. 

MOUNTED INSTRUMENTS 

Figure 10 shows computer plots of the instrument canister and wing 
mounts. The aft portion of the canister is simplified slightly compared 
with the detailed version shown in Figs. 5 and 6. The wing paneling is 
designed such that the mount plates (Fig. 10a) exactly cover two panels 
(i.e., the blackened panels in Fig. 9b). 

Figure 11 shows the instruments mounted under the Twin Otter wing. 

The instruments' axes make angles of 5° with the wing chord plane, forward 
end down. However, as mounted on the fuselage, the wing chord plane is 
tilted upward by 2.5° relative to the fuselage axis, so that the instruments' 
axes have a net downward tilt of 2.5° relative to the fuselage axis. 
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Figure 8. Digital description of the complete Twin Otter 

airplane suitable for use by the three-dimensional 
flow code. 
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b. View of underside 


Figure 9. Digital description of the Twin Otter used for this study. 

Panels covered by the instrument mount plates (Fig. 10) 
are blackened. Dimensional data and labeling are added to 
a computer plot. 
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b. PMS FSSP 


Figure 11. Computer plots of the wing-mounted instruments including 
short segments of the wing. Only the under side of the 
wing surfaces are plotted to avoid confusion of overlapping 
1 ines. 
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CALCULATION METHODS 


FLOW CALCULATION 

Air flow about the exterior of a three-dimensional body of arbitrary 
shape is calculated by a modified version of a first order panel code de- 
veloped by Hess (refs. 3, 9). The body surface must be approximated by 
contiguous, plane quadrilateral panels as illustrated by Figures 5-11 
above. Details of the method are given in references 3 and 9, and applica- 
tions are discuswed in references 1 and 2. 

With regard to accuracy, Hess (ref. 9) and Hess and Smith (ref. 10) 
report results of a large number of validation studies for flow about both 
lifting and nonlifting bodies. Panel code results are compared with exact 
results calculated from theory for simple bodies, and with experimental 
data for complicated bodies and combinations of bodies (e.g., wing-fuselage 
combinations). These comparisons show that the panel code accurately pre- 
dicts subsonic (i.e., essentially incompressible) exterior flows. On the 
other hand it is not capable of accurately predicting interior flows through 
tubes and ducts, and therefore cannot be used to predict air flux through 
the FSSP measurement tube. 

Air and hydrometeors must pass through the FSSP measurement tube, in- 
side of which the hydrometeors are sized by Mie scattering theory. The 
first order panel flow code discussed above can accommodate flow in/out of 
an orifice, providing flux in/out of the orifice is specified, but as noted 
it cannot accurately predict flow through the tube. A second order panel 
flow code is required to accurately simulate interior flows. Since a second 
order three-dimensional panel code was not readily available, some means 
was needed to estimate air flux through the FSSP measurement tube, which 
then could be input to the first order code. 

An axisymmetric , second order panel code (refs. 11,12) was used to 
calculate flow through the measurement tube for the measurement tube- 
canister combination of the isolated instrument; the measurement tube 
support arms (Fig. 6b) and airplane could not be included because of the 
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requirement that the body have axial symmetry. Average air flux into the 
orifice, which does not vary appreciably between 0° and 4° angles of attack, 
was calculated to be 95.9% of the free stream flux. To estimate the effect 
of the support arms on the flow, a three-dimensional calculation for the 

i 

canister-support arms combination, with the measurement tube omitted, 
and an axisymmetric calculation for the cannister alone were done. The 
difference in flow between these calculations at the intake orifice location 
j was 3.9% of the free stream flux, which again did not vary appreciably 

between 0° and 4° angles of attack. The net flux into the orifice, 

I 

accounting approximately for the effect of the support arms, is thus 
approximately 92% of the free stream flux. 

The three-dimensional first order code requires that the intake 
and exit orifices be paneled as though they are impervious surfaces and 
Figure 12 shows the paneling used: each orifice is represented by 

twenty-one panels. The flux through each orifice panel is specified in 
the input. For the FSSP in isolation, as well as when mounted under the 
Twin Otter wing with the airplane at 0° angle of attack, this specifica- 
i tion was taken to be 92% of the free stream flux as discussed above. 

However, calculation with the airplane at 4° angle of attack without the 
instrument shows an air flux at the location of the FSSP orifice of 
89% of the free stream value (owing to partial stagnation under the 
wing). Accordingly, with the airplane at 4° angle of attack, the flux 
through each panel was taken to be 82% (92 X .89) of the free stream 
flux. 


WATER DROP TRAJECTORY CALCULATION 

Assuming that the bulk air flow is not perturbed by the particles, 
and ignoring air density compared with water density, we use three dimen- 
sional, dimensionless equations of motion for water drops 


^ = 1 . 


dx N F 




D N R,s *"1 
D,s N R J 


N„N 
N 


( 1 ) 


I 
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Figure 12. Computer plots of the PMS FSSP measurement tube and 
its support arms. Intake and exit orifices of the 
measurement tube are each represented by 21 panels. 
(See Fig. 6 for dimensional data.) 
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Non-dimensional quantities are: 


V V a 


V Vr 


drop and air velocities 

still -air, terminal settling speed of the 
drop 

unit vector in the z (upward) direction 
time 

Davies number 


N F = V 2 /(Lg) 


Froude number 


N R = I ~ %| V Reynolds number 
C D water drop drag coefficient 


Dimensional quantities are: 


6 


equivalent sphere drop diameter 


p 


air density 


n 


air viscosity 


g gravity acceleration constant 

\I m free stream airspeed 

L a characteristic dimension of the body 


Here length is normalized by dividing by L, velocity by and time by 
L/V . N n and N D are for still-air, terminal drop settling (ref. 

13). 

We initiate the calculation far enough upstream to be essentially 
beyond the influence of the body where we can take v n = v = - £v c . We 
compute N R from these data, calculate Np from N R using experimental drag 
data for water drops (ref. 14), and proceed straightforwardly with a 
numerical integration of eq. (1) to obtain Vp and drop coordinates step- 
by-step from the free stream to the instrument. Air velocity, v g , is 
calculated via the flow codes as needed during the integration. Details 
are given in reference 3. 
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Accuracy of the numerical calculations has been studied in some de- 
tail (ref. 15), and trajectory results have been compared with experimental 
data where possible (refs. 3, 16). Results of these studies are summarized 
in reference 2. Comparisons with experimental data are satisfactory, and 
numerical integration errors are found to be negligible. 


CONCENTRATION FACTOR 


To estimate flow-induced drop flux distortions seen by the instruments, 
we compute a quantity called concentration factor, Cp, defined as 

r _ water drop flux at the instrument sampling volume /«\ 
l F water drop flux in the free stream ' ' 

Cp is estimated, for water drops of specified size, by computing a flux 
tube of drops from the unperturbed free stream to the instrument sampling 
volume. Figure 13 gives a crude illustration of such a flux tube. It is 
defined by a central trajectory (the dashed line in the figure) surrounded 
by a specified number* of trajectories that intersect the target plane 
approximately at the corners of a small, regular polygon about 
the center point of the instrument sampling volume. The initial plane is 
in the unperturbed free stream. These trajectories are computed by an 
iterative calculation that ensures that they intersect the target plane 
within a preset error tolerance of the desired target points.* 

Since the particle mass transfer rate through the tube is constant 
at all cross sections, it is easily shown that 


C 



(3) 


where A and A t are the flux tube cross sectional areas in the initial and 
target planes. 


*A sensitivity study (ref. 15) has shown that six flux tube boundary 
trajectories and a target point error tolerance of 25% of the target 
circle radius are quite adequate for concentration factor calculations. 
Accordingly we have used these parameter settings (NW = 6, TOL = 0.25; 
ref. 3, p. 84) for this study. 
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TARGET PLANE 



Figure 13. Artist's rendition of a water drop flux tube. The 
initial plane is normal to the central trajectory 
(dashed line), as is the target plane for the OAP. 

For the FSSP , the target plane is parallel to the 
measurement tube orifice plane. The initial plane is 
upstream in the unperturbed free stream. 
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For the OAP the flux tube center point in the target plane is at 
the center point of the laser beam (Fig. 5), and the flux tube boundary 
target points lie on a circle about the center point with a radius 
approximately equal to that of the laser beam. For the FSSP, the target 
plane is constrained to be parallel to that of the measurement tube 
orifice, and to lie a distance of 0.05" upstream of the orifice. Thus, 
the center of the flux tube in the target plane lies on the axis that 
passes through the center of the measurement tube, but lies a distance 

0. 05" upstream of the orifice center. 

For the FSSP calculations were done for two target circle radii: 

1. a radius slightly less (0.008" less) than the measurement tube inside 
radius (21/32"), and 2. a radius one half of the measurement tube inside 
radius. The coarser flux tubes are used to assess overall collection 
efficiency of the tube, while the finer flux tubes are used to assess 
collection efficiency for those droplets that are more likely to be accu- 
rately sensed and sized by the instrument. 

Trajectories of large drops are initiated at a distance of at least 
795" (10 fuselage diameters) upstream of the target planes. It is found 
that this distance is not always large enough for small droplets, par- 
ticularly for the wing-mounted FSSP at nonzero angle of attack, so this 
distance is doubled for droplets of diameters 20 ym.and less. 
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RESULTS 


FLIGHT CONDITIONS 

Since potential flow is calculated, the flow fields are independent 
of atmospheric properties and air speed, but trajectory calculations are 
dependent on these conditions. Usually the dependences are weak, but some 
dependence on air speed is found in this case, particularly for the FSSP. 

As expected, angle-of-attack is important for the wing-mounted instruments. 

All calculations are done for two angles of attack: 0° and 4° (forward 

end up) relative to the free stream direction (which direction is always 
normal to the gravity vector). In the case of an isolated instrument the 
angle is between the instrument axis and the free stream direction, while 
for an airplane-mounted instrument the angle is between the airplane fuse- 
lage axis and the free stream direction. Instrument axis angles of attack 
for the mounted instruments are discussed in the next section. 

Atmospheric properties are those at 7 kft in the U.S. Standard Atmos- 
phere (ref. 17). Specifically, air density is 0.993 kg m" and temperature 
is 274. 3°K. All calculations are done at two free stream air speeds: 95 

and 250 kts for the isolated instruments, and 95 and 130 kts for the 
airplane-mounted instruments. 

FLOW RESULTS 

Figure 14 shows the flow field for the PMS OAP mounted under the Twin 
Otter wing with the airplane at 0° angle of attack. Since the instrument 
axis makes a 5° angle (forward end down) with the wing chord plane, and 
the chord plane makes a 2.5° angle (leading edge up) with the fuselage 
axis, the instrument axis makes a net 2.5° angle (forward end down) with 
the free stream flow vector when the airplane is at 0° angle of attack. (For 
the 4° airplane angle of attack, the instrument axis makes a net 1.5° angle 
(forward end up) with the free stream flow vector.) Flow vectors at the 
center of the laser beam (shown with a circled root in Figure 14) have 
the following properties: 
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Figure 14. Flow field in the vertical plane normal to the OAP 
laser beam and through its center (•). The airplane 
fuselage axis is at 0° angle of attack to the free 
stream flow. 
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Airplane 
Angle of 
Attack 


Magnitude 
Relative to 
Free Stream 


Angle Relative 
to Horizontal 
Plane* 


Angle Relative 
to the Airplane 
Symmetry Plane + 


0 ° 


96.4% 


-9.2° 


1 . 1 ° 


4° 


83.1% 


-8.9° 


1 . 8 ° 


The center of the laser beam is approximately 10" from the wing surface. 

Figure 15 shows the flow field for the PMS FSSP mounted under the 
Twin Otter wing with the airplane at 0° angle of attack. Net instrument 
axis angles relative to the free stream vector for airplane angles of 
attack of 0° and 4° are the same as for the OAP. Figure 16 shows flow 
into the measurement tube orifice, again for an airplane angle of attack 
of 0°. Flow vectors at the center of the orifice have the following 
properties: 


Ai rplane 
Angle of 
Attack 


Magni tude 
Relative to 
Free Stream 


Angle Relative 
to Horizontal 
Plane* 


Angle Relative 
to the Airplane 
Symmetry P1ane + 


0 ° 


89.9% 


-12.7° 


2 . 0 ° 


4° 


76.7% 


- 9.7° 


3.1° 


The center of the orifice is approximately 11" from the wing surface. 


TRAJECTORY RESULTS 

We are interested in two quantities: 1. concentration factor (eqs. 

(2) and (3)) and 2. drop speed ratio, v p (eq. (1)). Deviation of concen- 
tration factor from unity expresses distortion in drop flux seen by the 
instrument. This distortion is caused by interaction of the drops with 
flow perturbations caused by passage of the airplane and instruments through 
the air. Instrument induced measurement biases are not considered here. 


★The horizontal plane is normal to the gravity vector, 

+These angles are such as to tilt the vectors to the outboard direction. 
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Figure 15. Flow field in a vertical plane parallel with the airplane 
fuselage axis and through the center point of the FSSP 
measurement tube intake orifice. The airplane fuselage 
axis is at 0° angle of attack to the free stream flow. 
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Figure 16. Flow into the FSSP measurement tube intake orifice. 

of the figure and angle of attack are the same as for 
Fig. 15. 



Drop speed ratio is the ratio of water drop speed to the free stream 
air speed* at either: the center point of the laser beam for the OAP, 

or into the center of the measurement tube orifice for the FSSP. This 
quantity is important for the OAP because it is needed to produce undistorted 
two-dimensional images of hydrometeors. While common practice is to assume 
that drop transit speed through the laser beam is the same as airplane air 
speed, we see below that actually transit speed varies with: instrument, 

drop diameter, angle of attack and air speed. 

In interpreting trajectory data it is important to understand the 

following generalities: 

1. Large drops, which have high inertia, tend to ignore local 
flow perturbations such that the large-drop tails of plots 
of concentration factor or speed ratio vs. drop diameter 
approach unity. 

2. Very small drops, which have small inertia, tend to closely 
follow the air flow such that the small-drop tails of plots 

of concentration factor vs. drop diameter show effects of flow 
divergence (Cp < 1) or flow convergence (i.e., streamline 
compression) (Cp > 1). Speed ratio asymptotically approaches 
that of the flow. 

3. Intermediate sized drops tend to show less predictable effects. 

For example, plots of concentration factor or speed ratio vs. 
drop diameter may show maxima or minima, sometimes sharply 
peaked, depending on geometry, angle of attack and sometimes 
air speed. 

Figures 17 through 21 show computer plots of portions of water drop 
flux tubes to the instruments. For the wing-mounted instruments all of 
the plots are for: 0° airplane angle of attack, 95 kts free stream air 


*For the wing-mounted instruments free stream speed corresponds to the 
airplane speed. 
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Figure 17. Flux tube of 20 ym diameter water drops to the isolated 
OAP with the instrument axis at 4° angle of attack to the 
free stream flow. V =130 kts. 

co 
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Figure 19. Flux tube of 100 mm diameter water drops to the wing- 
mounted OAP. Airplane angle of attack is 0°. 

V = 95 kts. 






Figure 21. Flux tube of 100 um diameter water drops to the 

wing-mounted FSSP. Airplane angle of attack is 0°. 
V = 95 kts. 
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speed and for water drops of TOO and 1000 urn diameters. One important 
consideration made clear by the plots is that trajectories to or near 
the instruments are in no danger of intersecting the wing surface. 

OAP 

Concentration factor results for the OAP are plotted versus water 
drop diameter in Figure 22, and Figure 23 shows a similar plot for speed 
ratio. Results for the isolated instrument and for the instrument mounted 
under the Twin Otter wing are included.* 

For the isolated instrument, calculations were done for free stream 
air speeds of 95 and 250 knots, and for 0° and 4° angles of attack for each 
speed. There is little variation of concentration factor and speed ratio 
with either air speed or angle of attack so that average values are plotted. 
The slight fall -off of concentration factor and speed ratio as drop size 
decreases is caused by the effect at the laser beam of flow stagnation 
against the forward hemisphere of the canister (Fig. 5). 

To fully understand results for the wing mounted instrument, it is 
necessary to more carefully consider flow conditions under the Twin Otter 
wing. Since at 0° angle of attack the wing is tilted (leading edge up) 

2.5° relative to the free stream, there is a slight stagnation condition, 
with attendant slight divergence and reduction in flow speed, under the 
wing in the region immediately forward of the instrument (Fig. 14). This 
condition primarily affects drops in the intermediate size range (diame- 
ters in the approximate range 20 to 200 pm). However, further aft this 
stagnation condition is overcome by slight convergent flow that is re- 
quired to negotiate the convex surface of the underside of the wing 
leading edge, and indeed, at the location of the laser beam, air flow 
speed is found to be slightly greater than free stream with the airplane 
at 0° angle of attack and without the instrument (Appendix A). With the 
airplane at 4° angle of attack, the under-wing partial stagnation condition 


*A11 concentration factor and speed ratio results are tabulated in 
Appendix A. 
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ISOLATED 



36 


Figure 22. Concentration factor (this study) or expansion ratio (Drummond) vs. water drop 
diameter for the PMS OAP. Average results for 0° and 4° angles of attack and 
V ro = 95 and 250 kts for the isolated instrument, and for V*. = 95 and 130 kts for 
the wing-mounted instruments are plotted. Drummond's results are taken from 



d31N3D WV38 d3SVl IV OllVd 033dS dOdO 
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Figure 23. Speed ratio vs. water drop diameter for the PMS OAP. Average results for 0° and 4 
angles of attack and = 95 and 250 kts for the isolated instrument, and for V ro 
= 95 and 130 kts for the wing-mounted instrument are plotted. Drummond's results 
are taken from Figs. 4 and 6 of ref. 6. 


is, of course, more pronounced and prevails both forward of the instrument 
and at the location of the laser beam. 

Flux and speed ratio results for the wing mounted OAP are completely 
consistent with this picture of the flow. For the airplane at 0° angle 
of attack the curves in Figs. 22 and 23 show minima at about 70 vm which 
are caused by the upstream stagnation effect on the midrange size drops. 

As drop size decreases, results rise slightly and level off slightly below 
unity principally because of the stagnation effect at the laser beam caused 
by divergent flow about the canister. For 4° angle of attack, the falloff 
of the curves with decreasing drop size is much more dramatic, as well as 
more persistent, as expected. These results predict undermeasurement by 
about 17% for both flux and drop speed throughout the cloud droplet size 
range, when the airplane is at 4° angle of attack. 

It is particularly interesting to compare Figs. 22 and 23 with Figs. 

24 and 25. (In Figs. 24 and 25 are plotted concentration factor and speed 
ratio vs. drop diameter at the location under the wing of the OAP laser 
beam, but with the instrument and its mount omitted.) This comparison 
indicates major contributions of the instrument itself to the distortions. 

At 0° angle of attack roughly half of the divergence effects for intermediate 
size drops, and all for small drops is caused by the instrument. At 4° 
angle of attack, slightly less than half of the effects for both inter- 
mediate and small drops are caused by the instrument. Implications of this 
are discussed below. 

Also plotted in Figures 22 and 23 are results of a similar study 
by Drummond (refs. 5, 6). The Canadian National Aeronautical Establish- 
ment has a Twin Otter with an OAP mounted under its wing at the same 
location as the NASA OAP (Fig. 9b), and the Drummond studies were done 
for this instrument. The only difference in the Canadian instrumentation 
is that the mounting pylon (Fig. 10) is 1" long compared with 5" for the 
NASA pylon. This puts the laser beam at a distance of approximately 
6" from the wing surface compared with 10" for the NASA instrument. 

Drummond used two-dimensional Joukowski airfoil theory to calculate 
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Figure 24. Concentration factor vs. water drop diameter at the location under the Twin 
Otter wing of the center of the OAP laser beam but without the instrument. 
Average results for = 95 and 130 kts are plotted. 


AOA 



40 


Figure 25. Speed ratio vs. water drop diameter at the location under the Twin Otter 
wing of the center of the OAP laser beam, but without the instrument. 
Average results for V = 95 and 130 kts are plotted. 


flow about the wing. Effects of the instrument on the flow were approxi- 
mated by adding the potential function contribution for a sphere to that 
of the airfoil, the sphere being placed at the location of the forward 
canister hemisphere (Figs. 5, 10 and 14), while the remainder of the in- 
strument structure was ignored. Expansion ratio (Fig. 22) is the two- 
dimensional analog of concentration factor. Drummond's water drop 
trajectory calculation methods are discussed in detail in reference 6. 

Drummond expresses flight conditions in terms of lift coefficient, 

C^. Using data and equations given by Drummond we convert to wing angle 
of attack to arrive at the following: 


C 


L 


Wing 

Angle of 
Attack 


Angle of Attack 
In Terms of 
Fuselage Axis 


0.348 


3.6° 


1 .1 


o 


0.784 


8 . 0 ° 


5.5 


o 


As already noted, when mounted on the fuselage, the wing has a 2.5° 

(leading edge up) angle of attack relative to the fuselage axis. In 
the rightmost column above we have subtracted 2.5° from the wing angle of 
attack so that it can be compared with angle of attack as defined here. 
Thus, Drummond's = 0.348 corresponds roughly with our 0° angle of 
attack, and his = 0.784 corresponds roughly with our 4° angle of attack. 
Comparisons of the corresponding curves in Figures 22 and 23 show good 
agreement between the results, in spite of the differences in geometry 
and calculation methods, except for the smallest droplets at the smaller 
angle of attack. For these cases Drummond's curves rise substantially 
higher than ours. This is probably caused by the fact that the Canadian 
instrument is 4" closer to the wing surface where the flow is likely to 
be substantially more convergent. 
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FSSP 


Results for the FSSP are plotted in Figs. 26-31. Concentration 
factors were calculated using flux tubes of two radii as described above 
on p. 23. Those calculated via the coarser flux tubes (solid curves in 
Figs. 26 and 27) predict overall collection efficiencies of the measure- 
ment tube, while those calculated via the finer flux tubes (broken curves 
in Figs. 26 and 27) predict collection efficiencies at the portion of the 
measurement volume inside of the measurement tube that is optically more 
active and that gives more accurate results. 

For the isolated instrument, results (Figs. 26-29) are insensitive 
to change of angle of attack, as was also found for the OAP, but con- 
trary to findings for the OAP, the results are slightly but sig- 
nificantly sensitive to change in air speed. For the instrument mounted 
under the Twin Otter wing, we find substantial sensitivity to angle of 
attack, as was also the case with the OAP, and in addition there is a 
consistent, though slight, sensitivity to air speed. Moreover, except 
for the largest drops, for which we expect little flux or speed distortion 
in any case, the distortions are more severe than for the OAP in all 
cases. Since differences in results at the two instrument locations 
without the instruments (compare Figs. 30 and 31 with Figs. 24 and 25) are 
minor*, it must follow that the bulk of the differences are caused by 
differences in flow around and through the instruments. This conclusion 
is particularly evident when results for the isolated instruments are 
compared. Thus, we find that the presence of the instruments, and even 
details of their geometries and operations, cause significant flow effects. 
This is an important finding since it has been common practice to assume 
that effects of flow perturbations caused by instruments themselves are not 
significant compared with effects caused by flow about the aircraft on 
which they are mounted. 


*The most pronounced differences are for the 0° angle of attack curves, 
which reflect differences in flow convergence under the wing at the 
measuring volumes of the two instruments. The OAP measuring volume is 
slightly closer to the wing surface, where there is slightly greater 
convergent flow, than is the case for the FSSP. 
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Figure 26. Concentration factor vs. water drop diameter for the PMS FSSP with Voo = 95 kts. 

Average results for 0° and 4° angles of attack are plotted for the isolated 
instrument. Flux tube radii at the measurement tube orifice are very slightly 
less than the orifice radius for the solid curves, and one half of the orifice 
radius for the broken curves. 
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Figure 28. Speed ratio vs. water drop diameter for the PMS FSSP with V*, = 95 kts. 

Average results for 0° and 4° angles of attack are plotted for the 
isolated instrument. 


I SOL ATE 0 V = 250 Kts 
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Figure 29. Speed ratio vs. water drop diameter for the PMS FSSP with \l m = 130 

and 250 kts respectively for the wing-mounted and isolated instrument. 
Average results for 0° and 4° angles of attack are plotted for the isolated 
instrument . 
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Figure 30. Concentration factor vs. water drop diameter at the location under the Twin 
Otter wing of the center point of the FSSP measurement tube intake orifice, 
but without the instrument. Average results for = 95 and 130 kts are 
plotted. 
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Figure 31. Speed ratio vs. water drop diameter at the location under the Twin Otter 
wing of the center point of the FSSP measurement tube intake orifice, but 
without the instrument. Average results for \l m = 95 and 130 kts are plotted. 


Regardless of whether we consider the solid or broken curves in Figs. 

26 and 27, we predict significant flux distortions for droplets in the 
cloud droplet size range (diameters < 50 ym). Quantitatively, these re- 
sults predict flux measurement errors by the Twin Otter FSSP on the order 
of 12 to 25% - caused by flow perturbations alone. The reasons for this 
are clear. Not only does the structure of the support arm-measurement 
tube combination present a considerable obstacle to the flow, but the 
air and drops must enter and traverse the measurement tube. 

As discussed above (pp. 17, 18) air flux through the measurement 
tube is reduced 8% relative to free stream for the isolated instrument, 
while for the wing mounted instrument, flux into the tube center is re- 
duced by 10 and 23% relative to free stream for 0° and 4° angles of attack 
(see the table on p. 26). The results in Figs. 26 through 29 are consistent 
with this. As noted previously, both concentration factor and speed ratio 
approach air flux (i.e., air speed) values as particle size approaches 
zero, whereas for large, massive particles both concentration factor and 
speed ratio are essentially unity. Thus, concentration factors calculated 
from the finer flux tubes (broken curves in Figs. 26 and 27) and speed 
ratios (Figs. 28 and 29) are expected to approach the values of 92% 

(isolated instrument), 90% (wing mounted, at 0° angle of attack) and 
77% (wing mounted, at 4° angle of attack) for the smallest droplets. 
(Concentration factor curves calculated from the coarser flux tubes are 
expected to approach lower minima for small droplets because the coarser 
tubes are more affected by diverging flow about the tube orifice edges 
and other obstructing structures.) These expectations are essentially 
realized for speed ratios in all cases, and for the concentration factors 
at 4° angle of attack for the wing mounted FSSP. On the other hand, 
concentration factor curves for the isolated and wing mounted FSSP at 
0° angle of attack pass through minima at about 7 ym droplet diameter and 
then rise fairly sharply toward the small-drop limit. It is surprising 
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for droplets as small as 2 ym that concentration factors remain considerably 
below the limits, though the data indicate that the limits would be finally 
reached by still smaller droplets. 

Preferred Orientation 


Preferred orientation of hydrometeors can be studied by examination 
of their drag vectors. It is well known that for Reynolds numbers large 
enough to indicate nonviscous flow but not so large as to indicate un- 
steady motion, hydrometeors orient such as to present maximum resistance 
to movement. Thus water drop flattening planes, ice column long axes and 
ice plate faces all orient normal to their drag vectors. Moreover, Beard 
(ref. 18) has shown that, at least for water drops, response time is short 
enough to allow reorientation to occur along curved trajectories at normal 
flight speeds. It is particularly important to know orientation of hydro- 
meteors as they pass through the measurement volume of a two-dimensional 
imaging OAP spectrometer so that the images can be properly interpreted. 
Water drop drag vectors at final points of flux tube central trajectories 
(i.e., at the centers of the instrument sampling volumes) as tabulated 
in Appendix B. These data predict that orientation varies significantly 
with all of: drop size, angle of attack and free stream air speed. More- 

over, the complexity of the data, particularly for the wing-mounted in- 
struments, indicates that additional study is warranted. 
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CONCLUSIONS 


ISOLATED INSTRUMENTS 

For the OAP we predict insensitivity of both concentration factor (i.e., 
flux ratio) and speed ratio* to both angle of attack and free stream air- 
speed. The only distortions are slight slowing and slight flux divergence 
for small droplets caused by an upstream effect, at the location of the 
laser beam, of flow stagnation against the forward hemispherical end of 
the canister (Fig. 5). These distortions are of little practical importance. 

The FSSP presents greater obstruction to flow than the OAP owing to the 
presence of the measurement tube and to the large size of its support arms 
(compare Figs. 5 and 6), and in addition requires that both air and drops 
traverse the interior of the measurement tube. Again no angle of attack 
effects are found, but persistent free stream air speed effects are found 
in the diameter range from about 5 to 50 pm for both concentration factor 
and drop speed ratio, though these effects are of no practical significance. 
Both concentration factor and speed ratio approach limits of about 0.9 
for small drops, which is consistent with our estimate of air flux through 
the measurement tube. The results predict that flux measurements of droplets 
in a size range typically found in clouds (diameters < 20 urn) would be too 
low by about 10%. Measurements of rain-)size drops are little affected by 
the flow. 

WING-MOUNTED INSTRUMENTS 

For the OAP, results are found to be insensitive to free stream air 
speed, but substantial angle of attack effects are found. At 0° angle of 
attack, distortions of drop flux and speed are small enough to be of 
little practical consequence. However, at 4° angle of attack sig- 
nificant decreases in both drop speed and flux are found for droplets with 
diameters less than about 150 ym. Much of this is cause by partial flow 


*Flux and speed ratios are for water drops at the instrument sampling volumes 
relative to free stream values. 
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stagnation against the under side of the uptilted wing. Minima in curves 
of concentration factor and speed ratio vs. drop diameter are 0.83 and 
0.81 respectively. The concentration factor curve remains constant at 
its minimum value for drop diameters less than 20 ym. The speed ratio 
curve reaches its minimum at about 30 ym, and then rises to stabilize at 
about 0.83 for smaller drops. These results predict undermeasurement by 
about 17% of both flux and transit speed through the laser beam for all 
cloud-size droplets, when the airplane is at 4° angle of attack. Our 
results agree well with the two-dimensional calculation results of Drummond 
(refs 5, 6). 

For the FSSP we again find very significant angle of attack effects 
on both drop flux and speed ratios and in all cases, the distortions are 
substantially greater than are found for the 0AP. Air speed effects are per- 
sistent for small drops, but are of no practical significance. For small 
droplets, in the cloud droplet size range, we predict undermeasurement of both 
flux and speed of about 10-13% at 0° angle of attack, to about 24% at 4° 
angle of attack. We repeat that these errors are caused by effects of 
flow perturbations alone. Measurements of rain-size drops are little 
affected by the flow 

Comparisons of results calculated with and without the instruments 
mounted under the Twin Otter wing show that roughly half of the flux and 
speed distortions are caused by the presence of the instruments themselves. 
This finding is quite significant because it has been common practice to 
ignore the presence of instruments on the assumption that flow perturbation 
effects caused by them are not likely to be signficant compared with those 
caused by the aircraft on which they are mounted. 

PREFERRED ORIENTATION 

Drop drag vector angles at the instrument measurement volumes appear 
to be rather complex functions of drop size, angle of attack and free 
stream air speed (Appendix B). This needs additional study because of 
its implications with regard to preferred orientations presented to the 
instruments of: flattened water drops, ice columns and ice plates. 
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APPENDIX A 


CONCENTRATION FACTOR AND SPEED RATIO DATA 
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APPENDIX B 


DRAG VECTOR ANGLES 

The drag vector, is given by the first term on the right side 
of eq. (1) (p. 18 ). Angle a is the angle between the projection of 5 in 
the z = 0 plane (i.e., horizontal plane) and the x axis vector, and g is 
the angle 15 makes with the z = 0 plane. The geometry is 


z(up) 



For the wing-mounted instruments we take the outboard direction (i.e., 
toward the wing tip) to be the direction of y. 

Tabulated drag vector angles (degrees) are for the final points of 
the central trajectories of flux tubes of water drops from the free stream 
to the instruments. For the OAP this point is close to the center of the 
laser beam, while for the FSSP it is close to the center of the measure- 
ment tube intake orifice. 
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